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Abstract 

If one introduces causality into quantum cosmology, then the prescription for the no-boundary 
universe should be revised. We show that the thermodynamic arrow of time associated with the 
perturbation modes should be reversed at the maximum expansion for the oscillating Hawking model. 
To an observer equipped with the time arrow, the universe will terminate its evolution after an half 
cycle. 
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In the no-boundary universe [1] , the Hawking model has been extensively investigated [2] . It is a 
closed FRW universe coupled to a scalar field <f> — <f)(t) with potential V{4>) — m 2 <f) 2 . Its Euclidean 
metric is described by 

ds 2 = d t 2 + b 2 (t) (d X 2 + sin 2 X (d0 2 + sin 2 6d^ 2 ) ) . ( 1 ) 

The Euclidean action is 

; = ^/^ + i-^-^). (2 , 

where dots denote derivatives with respect to the imaginary time t — it and the scalar field has 
been rescaled by multiplying the factor of (47r)2 for convenience. 
The field and the scale b obey the equations 

- 360 2 , 260 2 m 2 60 2 

+ —r~ = m <j>, b= — (3) 
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and the Hamiltonian constraint 



b 2 1 6 2 m 2 6 2 , s 

^-^-y + -r^°- (4) 

The quantum state of the universe is defined by the path integral over all compact 4-manifolds 
with the configuration of the wave function, the 3-geometry and matter field on it, as the only 
boundary. The wave function can be approximated by the exponential of the negative of the action 
of the complex solution with the boundary. This is called the WKB level. We shall work at this 
level in this paper. If one uses a Euclidean solution, then it can be described by a deformed 4-sphere. 
Its south hemisphere can be approximated by a half 4-sphere and the scalar field increases slowly 
from its south pole. As soon as the 4-sphere reaches its maximum size at the equator, the scalar 
field will increase rapidly and then the manifold will collapse into a singularity. 

For a general complex solution, the regularity conditions at the south pole ts implied by the 
no-boundary proposal are 

6 = 0, 6=1, = 0o, = 0. (5) 
The singular manifold is parametrized by the initial value 0o there. 
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A regular Euclidean solution is called an instanton. It has been shown that there does not 
exist any nontrivial regular instanton in this model [3] [4] . Equivalcntly, one cannot find a compact 
Euclidean regular solution with a 3-geometry (the equator) as the only boundary on which the 
second fundamental form vanishes and the normal derivative of the matter field is zero. At best, 
one can only find a Euclidean solution with approximately vanishing momenta at the equator [5]. 
Even when one relaxes the solution to be complex, the situation associated with the singularity will 
not change. 

One way out of the trouble caused by this singularity behavior in the scalar model is to reinterpret 
the Euclidean solution to the field equation as a constrained gravitational instanton [4] [6] . The 
south hemisphere of the manifold is the stationary action solution under the condition that the 
3-geometry is given, at the maximum size where the quantum transition is supposed to occur. The 
whole manifold is made by joining this south hemisphere and its oriented reversal as the north 
hemisphere. One can also use (f>o to parametrize the 3-geometry. The variational calculation shows 
that the stationary action solution should be regular and satisfy the field equations everywhere, with 
the only possible exception at the 3-geometry equator. Therefore the joint manifold has a stationary 
action under the restriction imposed at the equator or for a fixed 4>q and qualifies as a constrained 
gravitational instanton. The constrained instanton can be used as the seed for the creation of the 
universe. For this model, the canonical momentum of the scalar field at the equator, i.e, its normal 
derivative, is allowed to be nonzero. However, the nonzero real part of the momentum conjugate to 
the scalar in imaginary time is identified as its nonzero imaginary part in real time, therefore it can 
not lead to a real evolution in the Lorcntzian regime. 

The Lorentzian evolution of the model can be obtained through an analytical continuation at 
the equator. As mentioned above, in order to obtain a real evolution in the Lorentzian regime with 
real time one has to impose the condition at the equator that the imaginary parts of all fields and 
the real parts of all conjugate momenta in imaginary time are zero [7]. For the Hawking model this 
is 

Im(b) = 0, Im{4>) = 0, 
Re(b) = 0, Re{4>) = 0. (6) 
Apparently, to meet condition (6), one has to find a complex solution or complex constrained 
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instanton to replace the Euclidean south hemisphere. One has two degrees of freedom at the south 
pole, that is cf>Q e + i4>Q m , for the complex solution, and in addition to these one has two more degrees 
of freedom for the location of the equator in the complex time plane. At the mean time, there are 
4 conditions in Eq.(6) at the equator, therefore there are at most a discrete set of initial values 
of 4>q leading to Lorentzian evolutions. We set r\ = at the equator and then at the south pole 
Tg = —Tq c + jtq" 1 . Since the Lorentzian condition (6) at the equator is sufficient to guarantee the 
evolution along the line r Re — to be real, then at any point on this line, condition (6) is satisfied. 
One can further choose the point where Im(b) = as the equator. It means that at the creation the 
initial expansion rate is zero. This conclusion will not be changed even if one includes more matter 
fields into the model. 

To find the complex solutions leading to purely Lorentzian evolutions is interesting. But this is 
not the focus of this paper. Since there does not exist any regular instanton, at the equator the 
imaginary parts Im((f>) will not vanish. This means that in the Lorentzian regime the initial time 
derivative of the scalar at t\ is nonzero. However, one can select a suitable value 4>o such that Im{<p) 
is very small. After the creation the universe will undergo an inflationary period. Then the scalar 
4> decreases significantly and starts to oscillate. The inflationary period is succeeded by the matter- 
dominated phase of the big bang model. The universe will reach the maximum size at = — iT\ m 
and then recollapsc. 

In this paper we are particularly interested in the oscillating universe. When the universe reaches 
the maximum size at T2, it will meet condition (6) again. If at the maximum expansion the scalar 
field satisfies <p = or cj>' — 0, where prime denotes the derivative with respect to real time, then the 
contraction phase of the universe is exactly the time reversal of the expansion phase. For the first 
case the only difference is that the sign of the scalar should be changed. After the universe shrinks 
into the initial size, it will bounce back, and oscillate ad infinitum. The oscillating condition can 
be realized by adjusting the mass parameter m, in addition to the four parameters mentioned above. 

There are scalar and tensor perturbation modes around the minisuperspace background. The 
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tensor modes describes the primordial gravitational waves. They are 



2d„G" 



where G™- are the transverse traceless tensor harmonics and d n are the amplitudes. 
From the Wheelcr-DeWitt equation one can derive [8] 



d¥™ ] 1 
i- 



dt 6 3 



4"\ (8) 

where is the wave function for the mode. This is the Schrocdingcr equation for an oscillator 
with a time-dependent frequency v n w n/6. 
The classical evolution obeys 

< + — ± + (n 2 - IK = (9) 



or 



{bd n )" + ( n2 - 1 - y) M « = o, 



(10) 



where, and for the rest of the paper, the prime denotes the time derivative with respect to Lorentzian 
conformal time r\ defined by drj — b~ l dt. 

Now we try to work out the quantum state of the tensor mode using the no-boundary proposal. 
At the moment r = 0, one has the Lorentzian condition d n = and Im(a n ) = 0. These conditions 
can be satisfied by the choice of complex values of a n and d n at the south pole. The regularity 
condition at the south pole is b d n = 0. This condition can be realized by evolving Eq. (10) towards 
the south pole. Since one has b = (r — ts) + A(r — ts) + .--(A = const.) at the south pole, then 
the potential barrier is finite there, so bd n must be finite and b 2 d n is equal to zero at ts- In the 
Euclidean regime the term b" jb can be ignored if n 3> b. This implies that the complex solution is 
d n w d n (0) cosh(6 _1 nr). It is noted that the regularity condition of the modes at the south pole in 
Ref. [8] is too strict. At r = 0, taking account of the fact that the initial expansion rate is zero, one 
can obtain the wave function of the mode at the WKB level as in [8] 

4"Wp(-in& 2 4). (11) 
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That is, the no-boundary condition implies that each mode initially takes the minimum excitation 
state allowed by the Heisenberg uncertainty principle in quantum mechanics. 

However, when the wave length of the mode becomes equal to the horizon scale, then the adiabatic 
approximation breaks down. The wave function will then freeze. Quantum mechanically, the state is 
described by a squeezed vacuum, while classically the oscillator undergoes a parameter amplification. 
After the mode reenters the horizon in the matter dominated era, the adiabatic approximation 
(n 3> b) will be recovered. However, the oscillator then is in a highly excited mode. It was argued 
that at the expansion phase the entropy associated with the perturbation modes increases, and the 
thermodynamic arrow of time is implied by the cosmological arrow which is defined by the direction 
of the expansion of the universe [9] . 

Since the operator in the right hand side of Eq. (8) is time symmetric, the amplitude of the 
wave function with initial value (11) is an even function of time at the creation moment. For the 
oscillating model, if we identify the two ends of the Lorentzian time in a cycle, the amplitude of 
the wave function will remain symmetric about the identified time origin. The amplitude evolved 
forward from the big bang and that evolved backward from the big crunch will be equal at the 
maximum expansion. It is noted the Schroedinger equation is of first order in time derivative. If 
one follows the evolution from the big bang until the big crunch, he will find that the wave function 
will return to the initial state at the crunch end up to a phase. 

The classical equation (9) can be identified as that for a particle tunneling through a potential 
barrier 1 — n 2 + b" jb in one dimensional quantum mechanics with coordinate r\. In general, the 
solutions are not symmetric about the maximum expansion due to reflection of the potential barrier. 

(n) 

However, the no-boundary wave function tyj, should be the quantum counterpart of a classical orbit 
which is symmetric about the maximum expansion. The time arrow is determined by the direction 
of increasing entropy and the highly excited perturbation modes are of higher entropy, therefore the 
arrow of time should be reversed at the maximum expansion. 
The scalar perturbation of the matter is 




n 



(12) 
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where Q n are scalar harmonics. The scalar perturbation of the metric is 



E 



\/6 I fc„p^ 2a n fl ij Q n + 66„P£ 



(13) 



where P™ = Q™ n /(n 2 — 1) and P/j = f^-<3™/3 + Q^j/in 2 — 1). fc„ and <?„ are Lagrange multipliers 
and induce two constraints. Therefore, there is only one true scalar degree of freedom. In the 
b n = k n = 0,g n = — a n gauge, one has the constraint [10] 



a' n + ja n = -ty'fn 



and the decoupled equation of motion for a n [10] 
a" + 2 



2ll 



2b', 



+ (n 2 + 3) 



a n = 0, 



(14) 



(15) 



which can be rewritten as the equation for a n = ba n /(p' 



b" 

~ + 2| --r 



in) ^' + (- 2 + 3) 



o. 



(16) 



Again, this takes the form for a particle tunneling through a potential barrier in one-dimensional 
quantum mechanics. The quantity a n behaves like the wave function. Form (16) is very useful, 
since the operator within the square bracket is of even order in the time derivative, and the dissi- 
pation effect disappears. It can be continued into the Euclidean regime by the Wick rotation. The 
Lorentzian condition for a n at the equator is a n — and Ima n = 0. For the case n^> b, all terms in 
the square bracket of Eq. (16) are negligible in comparison with (n 2 + 3) in the Euclidean regime. 
One can use a similar method to obtain the complex solution to meet the Lorentzian condition at 
the equator and the regularity condition at the south pole as in the case of the tensor modes. 

The above solutions can be used for construction of the wave function for the scalar mode at the 
creation moment and obtaining the minimum excitation state similar to Eq. (11). For simplicity, one 
uses the oscillating model with (f> = at the maximum expansion. By identifying the big bang and 
the big crunch ends Eq. (16) becomes symmetric about the creation moment. The Schroedinger 
equation corresponding to its classical counterpart (16) can be written. Then we use the same 
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argument as for the tensor modes to get a wave function with time symmetric amplitude at the 
minimum size and the maximum expansion. Therefore, the time arrow associated with the scalar 
mode also flips at the maximum expansion. 

From the no-boundary philosophy, the wave function for the perturbation modes is a function 
of the configuration only and should be independent of the orbit of the minisupcrspace background. 
In this way, the quantum state of the modes at the big bang end should be identical to that at the 
big crunch end. However, when one uses analytic continuation at the WKB level, one may worry 
that the quantum state of the modes may depend on the contour. However, it turns out that this is 
not the case. 

In summary, for both the tensor and scalar cases, the perturbation modes start as a minimum 
excitation at the big bang end, it will return to the same state at the crunch end. The entropy will 
increase when the modes evolve from the minimum excitation to the highly excited state. Thus, the 
thermodynamic and psychological arrows of time should be flipped at the maximum expansion [9]. 

It is noted that the appearance of singular terms encountered in Eq. (16) is superficial. It has no 
effect in Euclidean regime. In the Lorcntzian regime, all these singular behaviors must be cancelled 
or one can simply integrate the equation in the complex time plane along a very close time line 
parallel to the Lorentzian contour and then take a limit to approach the true Lorentzian evolution. 

Our conclusion about the scalar modes is different from Rcf. [10]. The main reason is that the 
classical Lorentzian evolution can be described by the wave function of a particle penetrating through 
the potential barrier. If one lets the wave function at the big bang end be a pure propagating wave 
as in [10], then at the other end, it must be a superposition of propagating and reflecting waves. 
There is no surprise that the amplitude at the crunch end becomes much larger! 

To set the pure propagating wave at the big bang end is wrong. First, the complex values of d n 
and a n are unphysical. Second, and more importantly, for the whole perturbation calculation leading 
to (8) one has not taken account of causality condition. The naive analytic continuation leads to 
complex d n and a n . In quantum field theory one has to decompose the field into the positive and 
negative frequency parts associated with creation and annihilation operators [11]. The Lorentzian 
condition complies with the causality condition. This requirement has no effect on the background 
calculation, since one ignores causality in the minisupcrspace calculation anyway. 
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To recover causality, one has to take the real part of the whole wave function Re(C expiS), then 
one obtains a couple of complex conjugate Schroedinger equations for perturbation modes, they are 
for negative and positive frequency parts, respectively. The Hermitian property of the perturbation 
field leads to real values of d n and a n . This phenomenon was previously encountered as dealing 
with the perturbation modes around the Schwarzschild-de Sitter background [12]. The modified 
prescription will leave the prediction on the origin of the structure in the universe intact, but it does 
change the issue about the time arrow in the universe. Indeed, if one chooses the pure propagating 
wave at the big bang end, one has implicitly introduced the time arrow by hand from the beginning. 

The universe may keep oscillating after the big crunch. However, to an observer equipped with 
the thermodynamic or psychological arrow of time, he will find that the evolution will terminate at 
the maximum expansion. The contraction phase of the universe is identical to its expansion phase. 
One ends with a half-cycling universe [13]. 

For non-cycling Lorentzian evolutions, the above argument is no longer valid. As we know, one 
has to fine tune the parameters to yield the oscillating models. However, it is believed that the evolu- 
tions of the background and the quantum state of the modes are not so sensitive to these parameters 
before the anti-inflationary period near the big crunch. Therefore, for the general Lorentzian orbits 
the time arrow should be reversed at the maximum expansion as well, it might even be flipped again 
near the big crunch if the universe recollapses into a true singularity. However, even for an eternal 
observer, he would not be able to cross the maximum expansion to experience this. It is noted that, 
since the entropy reaches its maximum at the maximum expansion, the time arrow becomes very 
vague then. 

One can even present a conjecture: under the regularity condition at the south pole and the 
Lorentzian condition at the equator, all Lorentzian orbits of the background are oscillating or quasi- 
oscillating. Then what we argued for the oscillating orbits about the time arrow applies to all orbits. 
But I cannot offer a proof. 
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